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Analytical Study of Gas 
Turbine Combustor 


Abstract: The primary purpose of this study is to discuss the main 
requirements of gas turbine combustors in industrial engines. The 
principal geometric and aerodynamic features that are common to 
most types of gas turbine combustors are briefly reviewed. Reference 
is made to most of the key issues involved in combustor design and 
development. In this project, analytical formulae described by various 
researchers are studied and also their parametric study on the design 
of combustor is discussed in this report. The various combustor 
components such as diffuser, swirler, combustor liner, combustor 
casing and different zones are systematically designed and c- 
programming is provided for quick access of the design parameters. 
Manual calculations are somewhat tedious, so in order to save time 
and money this attempt is made. A technical discussion on combustor 
design procedure by using C programming is provided. 


Index terms: Combustor, Reference Area, Diffuser, Swirler, Liner, Liner 


holes. 


I. INTRODUCTION 

A combustor is a part of a gas turbine engine where 
combustion takes place. It is also known as a burner, 
combustion chamber or flame holder. In a gas turbine 
engine, the combustor is fed with high pressure air by the 
compression system i.e. compressor. The combustor then 
heats this air at constant pressure. After heating, air passes 
from the combustor to the turbine. 

The engines used a single chamber known as a can type 
combustor. Today three main configurations exist: can, 
annular and cannular. 


II. COMBUSTOR TERMINOLOGY 

The case is the outer shell of the combustor. The 
casing serves as a pressure vessel that must withstand the 
difference between the high pressures inside the combustor 
and the lower pressure outside. Diffuser is provided to slow 
the high speed, highly compressed air from the compressor 
to a velocity optimal for the combustor. The liner contains 
the combustion process and introduces the various air flow 
into the combustion zone The swirler are the part of the 
combustor that the primary air flows through as it enters 
the combustion zone. 

Their role is to generate turbulence in the flow to 
rapidly mix the air with fuel. The function of primary zone 
is to sustain flame & to rovide optimum temperature. Gas 
diluted with air through intermediate holes. This air 
induces further turbulence and supplies more oxygen to 
complete the burning of soot formed during primary 
combustion. Gas is further diluted by air admission through 
the dilution holes to reduce temperature to that level which 
is acceptable to the turbine blades. 
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Figure. 1: Combustor components 


In axial-flow compressors, the stage pressure rise is 
very dependent on the axial flow velocity. To achieve the 
design pressure ratio in the minimum number of stages, a 
high axial velocity is essential, in many gas turbine 
engines, compressor outlet velocities may reach 150 m/s or 
higher. It is impractical to attempt to burn fuels in air 
flowing at such high velocities. Here the fundamental 
pressure loss would be excessive. Thus, before combustion 
can proceed, the air velocity must be greatly reduced, 
usually to about one-fifth of the compressor outlet velocity. 
This reduction in velocity is accomplished by fitting a 
diffuser between the compressor outlet and the combustor 
inlet. 

a. Diffuser Geometry: 

The geometry of straight-walled diffusers may be 
defined in terms of three geometric parameters, as shown 
in Figure 2. 

1. Area ratio(AR) 
2. Characteristic length to inlet dimension ratio(L/R1) 
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3. Divergence angle(@) 
These are related to each other by relation 


= dee | oi ay2 
AR = 1+2-sin®@ + (7 sin®) .. (1.1) 


yt 
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Figure2: Diffuser Geometry 





1.Area ratio(AR): Function is to achieve prescribed 
reduction in velocity. 

AR=A2/Al1 
Continuity eq" ; A2xV2 = AlxV1 

A2/A1=V1/V2 
i.e. there is inverse variation in area ratio & velocity ratio. 
For max. AR there is more velocity reduction in diffuser 
i.e. high pressure recovery. 
2. Non dimensional length (L/R1): In combination with the 
area ratio the length defines the overall pressure gradient 
which is the principal factor in boundary-layer 
development. Usually, either the wall length, L, is used as a 
characteristic length. It is expressed in non dimensional 
form by dividing by a representative inlet dimension. 
3. Divergence angle(8) : In long diffusers of low 0, the 
pressure loss is high. With an increase in 0, both diffuser 
length and friction losses are reduced, but stall losses 
become more significant. 


b. Swirler Design: 

The primary-zone airflow pattern is of prime 
importance to flame stability. Many different types of 
airflow patterns are employed, but one feature common to 
all is the creation of a toroidal flow reversal that 
recirculates a portion of the hot combustion products to 
mix with the incoming air and fuel. These vortices are 
continually refreshed by air admitted through holes pierced 
in the liner walls, supplemented in most cases by air 
flowing through 
atomization. 


swirlers and by air employed in 

One of the most effective ways of inducing flow 
recirculation in the primary zone is to fit a swirler in the 
dome around the fuel injector. Vortex breakdown is a well- 
known phenomenon in swirling flows; it causes 
recirculation in the core region when the amount of rotation 
imparted to the flow is high. This type of recirculation 


provides better mixing than is normally obtained by other 
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means, such as bluff bodies, because swirl components 
produce strong shear regions, high turbulence, and rapid 
mixing rates. These characteristics of swirling flows have 
long been recognized and have been used in many practical 
combustion devices to control the stability and intensity of 
combustion and the size and shape of the flame region. 





Figure 3: Axial swirler 

The conventional notation for axial swirlers is indicated 
in Fig 3.2. This figure shows a flat-vaned swirler whose 
vane angle is constant and equal to 0, With curved-vane 
swirlers, the inlet blade angle is zero and the outlet angle is 
0,. 

An important design requirement is that the swirler 
should pass the desired airflow rate for a given pressure 
drop APsw, which is usually assumed to be equal to the 
liner pressure drop, AP;. We have , 


0.5 
2p3APsw 


m., = 
sw Kaw l 


(1.2) 





sec 0 1 
_ 


Asw AL 





(sec 6)2 
2p34Psw fs 1 
msw “Ksw Ay2 


— 2 Osw —Dnub* 
3 (Dsw —Dhub )? 
Thus, for a simple axial swirler, the minimum vane 
angle required to obtain strong recirculation (SN > 0.6) for 
a typical swirler having Dhub/Dsw = 0.5 is calculated from 
Equation 1.4 as 38°. 





(1.3) 


Sn tanO, ...... (1.4) 


c. Combustor Liner 

Combustor casing is provided with a liner. The 
combustor liner contains the combustion process and 
introduces the various airflows into the combustion zone. 
To establish air flow pattern within the liner, there is need 
to determine the number, size, shape and deposition of liner 
holes. Necessary for easy light-up, efficient and stable 
combustion, adequate wall cooling and delivery of hot 
gases with a suitable temperature profile. 

If liner wall contains a row of n holes, each of which 
has an effective diameter dj , then the total mass flow rate 
m, of air through these holes is given by 


my Som? aU asctesdeciteccases (1.5) 


14 
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Now 
0.5 
_— (2APL 
U; = (=) edna taeda te eeu ses (1.6) 
Hence, 
m, = an Gi(Zp30P >. Janie (1.7) 
Add, 
15.25 
1 P lececctces (1.8) 
(P3 Llp ) 
d; 
d, = e Res Vavsnusaraenraeses (1.9) 
D 


Ill. ALGORITHMS 


T3, P3, Ts, P;,Ts, AIF, C,, Testis Coe, R, e, f, ’ B, Py, A, Ca, 
h, th, Kyw these are floating type global variables. 


a. Algorithm to Calculate Mass Flow Rates and 
Reference Area 

This algorithm is used to calculate mass flow of fuel, mass 

flow rate of air and reference area of combustor. 


Pdef, ma3, ma, mf, mg, mpz, T;, Ts, W, W,, Pr, a, b, s, J, 
L, Cc, qref, arf, ar, Asefs Urers R,, Ay D, Las Li, Ly PF, 
Mopzms Mpzhs Maz, Miz Uhpz, Aipzs X, Zi, J, Msy are the various 
float type variables. 
START 
Read (Pressure ratio Pr) 
Ty 
pe ie 
Ts < Ty X Tots X (Ty — Ts) 
Write (Turbine exit temperature, Ts, Ts ) 
Read (Total work output in kW, W) 
W, = 2x W 


Ts <— 


, W, x 102 
mM! e¢ errr 
[(afr + 1) x C,, x (T, — Ts )] 


m, <— Mm, X afr 
M, 
<_ 
7 0.5 
e Mpz 
3 0.25 


b R 
<e_ 
2 


2 
m,, X 73 
c< | ——— 
Ps 
Write(Work required to run the turbine ,W,) 
Write(Mass flow rates, ms,m,,m,3) 


Aver — V(b X eX f X c) 


Ree 
a 
° RXxTs 


m 
Over 


Mp 





mM, 


a3 
(ee 

Ro x Aref 
Write(Reference area, A;er) 
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Write(Reference velocity, User) 
Read(Area ratio, ar) 
A; — ar X Aver 

A 
(P/4) 
Ly, <— 1.5 x Dy 
L,, — 0.5 x D, 
Paer — P3 — Py 
Py — Piep X 1075 
AP, = B X Puer 
dpl — AP, x 107° 


(3.12) 





Di <— 


(3.14) 


AP, 

Jeexfx— 
Ps 

Jie J x 107° 
Write(Liner area, A)) 
Write(Liner diameter, D)) 
Write(The liner lengths are Lg,, Liz) 
Write(The value of pressure differences, Pg, dpl) 
Write(The value of pressure drop, Ji) 


Read(The value of pattern factor, PF) 


wi 
<_ ud 
ONT = PE 





AXJxXs 

Loz = L on (Laz + Liz) 

Write(The liner lengths are, L;, Lp,) 

HALT 

functionTwo: 

b. Algorithm to Calculate Number of Holes on Liner 
This algorithm is used to calculate diameter of liner holes 
and number of liner holes. 

roundnpz, roundniz, roundndz are the various intiger 
variables. 

M,3, Ma, D), delpl, Ma, Mpzm, Msws Mpzhs Maz, Miz, Anpz, ipz, & 
Npz» Nizs Ajizs Anizs Anaz, diaz, Naz are the various float variables. 
START 

Mg — M3 X 0.75 

m 


pam “— My 


My <— mM 


pz 
Mg, — M,, X 0.3 


pzm 


mj, — m,, X 0.15 

Write (The mass flow rates are, Mg, Mpzms Mpzh, Maz, Miz) 
dnpz — 0.1 X Dy 

dn dipe VCd 

g <— yP3 x (delpl/T3) 

15.25. X Mpph 


ie 2 
diz xX 2 


pz — 


roundnpz < round(n,,) 
niz — Ty f 2 


roundniz < round(n,,) 
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15.25 x m;, 
dz. |S 
ni, X g 


di 
vCd 
dna, < 0.15 x D, 


diaz - dnaz x vVCd 
15.25 x mg, 
dian xg 





hin — 


Ng, — 


roundndz < round(ng,) 

Write(The hole diameters, dip, diizs diaz) 
Write(The number of holes, ny, niz, Naz) 
HALT 


c. Algorithm to round of the number of hole value. 
START 

tempRound is an integer value. 
calculatedDiffrence is a float value. 

int tempRound < (int)n,, 
calculatedDiffrence — n,, — tempRound 
if(calculatedDiffrence>=0.5) 

return tempRound+1 

else 

return 

tempRound 

HALT 

Function Three 


d. Algorithm to Calculate Swirler Diamentions 

This algorithm is used to calculate swirl number, swirler 
area and swirler diameter. 

Ry, m,, Aj, Delpsw, delpsw, i, u, Sy, j, Asw, Ny, mM, 0, g, Dyy, 
X_vane, x these are some float type variables. 

START 

Delpsw < P; x 0.03 

delpsw < Delpsw x Ore 

Write(The pressure difference across the swirler, delpsw) 
Read(The value of vane angle, x_vane) 

i — tan(0.0174 x x_vane) 





Dhub : 
S = cs ~— tan(0 
ae eo a 
_ 2XR, X Delpsw 
LS 

K,,, Xm,” 

1 + tan?(0v 

oe (0v) 






2X Ry X AP ay + i 
K,,, X m,2 A; 
Write(Swirl number, Sn) 
Write(The swirler area, Asw) 
Read (The number of vanes, nv) 
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m+« —(0.5 Xn, Xt, X (1 —h)) 
0 < (pi/4) x (1 —h’) 
Qg— m? a (4 xO X (—A,,)) 


—m+./q 


Dy = 

: 2Xo 
Write(The swirler diameter, D,,) 
HALT 


Function Four 


e. Algorithm to Calculate Diffuser Dimensions 

This algorithm is used to calculate diffuser inlet area, 
diameter, radius and diffuser length. 

Aj, Ao, V1, V2, AR, Ro, ma, Ay, Di, Ry, deffusion_angle, L, 
M, y, Ares, Ao=Ares, deffusion_angle=y these are some float 
variables. 

START 

Read(The diffuser inlet velocity, v,) 

Read(Diffuser exit velocity, v2) 


Vi 
AR<« — 
V2 
AAs 
é a 
1 2 Vv 
Ay 
D ae! 
' pi/4 
R Pi 
see 
SD 


Write(The value of area ratio, AR) 
Write(The diffuser inlet area, A;) 
Write(The diffuser inlet diameter, D,) 
Write(The diffuser inlet radius, R;) 
Read(The diffusion angle in degree, y) 
M < sin(0.0174 x y) 
—1+¥VAR 

* M/R; 
Write(Length of diffuser, L) 
HALT 
Function Five 


f. Algorithm to Calculate Fuel Injector Parameters 
This algorithm is used to calculate coefficient of discharge, 
flow number and orifice diameter. 


Delpf, Cd_inj, Do, Rs t, FN, angle, x_spr, ld, Z, m; these 
are some float variables. 

START 

Read(The value of length to diameter ratio, Id) 
Read(The fuel density, R,) 

Read(The fuel pressure difference, Delpf) 

Cdinj — 0.827 — 0.0085 x Id 


Ze IR, x Delpf x 10° 


ENS 
<_ — 
Z 
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1.11 x Cd_inj x Z 





Write(The coefficient of discharge, Cd_inj) 

Write(The flow number, FN) 

Write(The orifice diameter, D,) 

Read(The value of spray cone angle, x_spr) 
0.00805 x ./R; x FN 

< D, X cos(x_spr X 0.0174) 
Write(The film thickness, t) 
HALT. 


IV. RESULTS 
1. Effect of diffusion angle on diffuser length: 


(m) 
——) 
KR WN 


fmt 
vo 
e 
3 
= 


10 20 
Diffusion Angle(60) 





Figure 4: Effect of diffusion angle(0) on diffuser 
length(L_diff) 

As the diffusion angle increases the length of diffuser 
get reduce, which is beneficial in order to reduce the 
weight of engine. Again large diffuser deflection is 
required for large recirculation area and to reduce the 
pressure losses. 


2. Effect of various vane angle on swirl number: 





Swirl No.(Sn) 


50 





Vane angle(Ov 0) 


Figure 5: Effect of various vane angle(0,) on swirl 
number(Sn). 

As the swirl angle increases the swirl number also 
increases. For strong recirculation the swirl number should 
be greater than the value .0.6(Sn>0.6), for that proper swirl 
angle selection is important, after calculations it is clear 
that the swirl angle should be greater than 35° to get Sn>0.6 
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3. Effect of turbine work output on fuel mass flow rate 
distribution for different pressure ratios. 





70 
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Figure 6: Fuel mass flow rate(mf) Vs turbine work 
output(Wt) for different pressure ratios. 


As we increase the pressure ratio, the fuel consumption for 
different powers also increases 


4. Combustor Liner Length and Holes Distribution: 


oo 


Liner zone lengths(m) 
CoOorRNwW kN DW ~I 





Wwt(Mw) 
Figure 7: Liner Length (L) distribution Vs Turbine 
Power(W,) 

Length of liner for primary zone is larger than other 
two, as primary combustion takes place in this zone and air 
requirement is more for combustion. Length of dilution 
zone is next shorter where cooling of hot gases takes place 
to make them sustainable by turbine blades and next is 
intermediate zone where combustion incomplete 
combustion product takes place. 





Number of holes(N) 


Turbine power output(Wt) in Mw 
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Figure 8: Number of holes(N) distribution Vs turbine 

power output(Wt) at pressure ratio=3. 

For good cooling effect we need to increase the number 
of holes on the dilution zone, it is possible by reducing the 
diameter of dilution holes, because temperature at the 
combustor exit should required to be less so that turbine 
blades are able to sustain. 


5. Fuel Injector Parameters 
45 

40 

35 
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20 Pr=5 
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10 
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Figure 9 : Orifice diameter(D,) Vs turbine power 
output(Wt) for various pressure ratios. 
The nozzle diameter reduces with increase in pressure 
ratio. 


Fuel film thickness, 
tx10-2 (cm) 
OoOorRFNwW KN WANA CO OO 





Figure 10: Fuel film thickness(t) Vs turbine power 
output(Wt) for various pressure ratios. 

The decrease in fuel film thickness takes place as an 
increment in the pressure as well as turbine power output is 
done. 

6. Adiabatic Flame Temperature 





Adiabatic flame 
temperature(Tabd) 


0.5 1 


Equivalence ratio(@) 





Figure 11 : Adiabatic flame temperature(Tabd) Vs 
equivalence ratio(@). 
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At equivalent ratio is one(@=1), i.e. stoichometric 
composition, the adiabatic flame temperature reaches the 
peak temperature e.g. 2573K for Kerosene fuel. As the 
equivalence ratio increases, the adiabatic flame temperature 
also get increase upto stoichometric composition after that 
the temperature get reduces. 


IV. CONCLUSION 

As the manual calculations are very tedious and time 
consuming job. It should get overcome with the help of 
computer simulations for design of combustor. It is 
advantageous as in this practice no need to perform the 
experiment again and again user get the required 
parameters affecting combustor design by just providing 
required output from the gas turbine in the c program. It is 
also beneficial to select design giving better performance 
due to ease in comparison between two or more design. 


REFERENCES 

1. Lefebvre A.H.( 1965), “Progress And Problems In Gas- 
turbine Combustion” , Tenth Symposium (International) on 
Combustion, pp. 1129-1137. 

2. Gupta A.K.(1997), “Gas Turbine Combustion : Prospects 
And Challenges” ,Energy Convers, Vol. 38, No. 10- 
13,pp.1311-1318 . 

3. Gupta, A. K., Lilley, D. G., and Syred, N. (1984), “Swirl 
Flows”, Abacus Press, Tunbridge Wells, UK. 

4. Sovran, G., and Klomp, E. D.( 1967), “Experimentally 
Determined Optimum Geometries for Rectilinear Diffusers 
with Rectangular, Conical or Annular Cross Section, in G. 
Sovran, ed., Fluid Mechanics of Internal Flow, Elsevier, 
Amsterdam, pp. 270-319. 

5. Channwala S. A. & Kulshreshtha D.B. (2006),““Numerical 
Simulation to study the effect of divergence angle on 
performance of inlet diffuser section of gas turbine 
chamber’, Tenth International Congress of Fluid Dynamics 
,December 16-19, Stella Di Mare Sea Club Hotel, Ain 
Soukhna, Red Sea, Egypt, ICFD10-EG-3013. 

Books: 

1. Lefebvre, A. H. & Ballal, D. R.,2010. Gas Turbine 
Combustion Alternative Fuels and Emissions. 3" ed. CRC 
Press is an imprint of the Taylor & Francis Group,,. 

2. Samanta, D.,2006. Classic Data Structure. New Delhi, 7 
printing 2006. Prentice Hall of India Private Limited. 

3. Rajaraman, V., 1994. Computer Programming in C. New 
Delhi, on Edition, Prentice Hall of India Private Limited. 

4. Ganesan, V.,2000. Computer Simulation Of Compression- 
Ignition Engine Processes. Hyderabad, First published 
2000, Universities Press(India) Limited. 

5. Cohen, H., Rogers, GEF.C. &  Saravanamuttoo, 
H.LH.,1987. Gas Turbine Theory. 3"! Edition, Longman 
scientific and Technical , New York, 40-200. 

6. Mathur, M. L., Sharma, R. P., 1976. Internal Combustion 
Engine. 1" Edition, 23" Reprint 2012, Dhanpat Rai 
Publications (p) LTD, 894-923. 


www.rame.org.in 


